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Abstract 

Recently, Snider et al (2020 Nature 586 373) reported on the observation of 
superconductivity in highly-compressed carbonaceous sulfur hydride, Hx(S,C)y. The highest 
critical temperature in Hx(S,C)y by 5 K exceeds previous record of T. = 280 K reported by 
Somayazulu et al (2019 Phys. Rev. Lett. 122 027001) for highly-compressed LaHıo. In this 
paper we analyse experimental temperature dependent magnetoresistance data, R(T,B), 
reported by Snider et al. The analysis shows that Hx(S,C)y compound exhibited Te = 190 K (P 
= 210 GPa) has the electron-phonon coupling constant Ae-pn = 2.0 and the ratio of critical 
temperature, Te, to the Fermi temperature, Tr, in the range of 0.011 < T./Tr < 0.018. These 
deduced values are very close to ones reported for H3S at P = 155-165 GPa (Drozdov et al 
2015 Nature 525 73). This means that in all considered scenarios the carbonaceous sulfur 
hydride 190 K superconductor falls into unconventional superconductors band in the Uemura 
plot, where all other highly-compressed super-hydride/deuterides are located. It should be 
noted that that our analysis shows that all raw R(7,B) datasets for Hx(S,C)y samples for which 
Snider et al (2020 Nature 586 373) reported Te > 200 K cannot be characterised as reliable 
data sources. Thus, independent experimental confirmation/disprove high-T; values in the 


carbonaceous sulfur hydride is required. 


I. Introduction 

In 2015 Drozdov et al [1] discovered the first near-room-temperature (NRT) 
superconductor highly-compressed sulphur hydride, H3S. To date, NRT superconductivity 
has been observed in four super-hydrides/deuterides systems subjected to high pressure: Th- 
H [2], S-(H,D) [1,3-5], Y-H [6,7], La-(H,D) [8-10]. 

Recently, Snider et al [11] reported on the observation that ternary compound H,(S,C)y 
exhibits the transition temperature within a range of T. = 275-287 K (at P = 258-270 GPa), 
which is by about 5 K higher than previous record of T. = 280 K, reported by Somayazulu et 
al [8] in highly-compressed LaHio9. While detailed phase/structural and phonon spectrum and 
other physical properties measurements [1-7,9,12-14], as well as the first-principle 
calculation studies [2,6,15-27] are on-going tasks for Hx(S,C)y compound [28], in this paper 
we report results of the analysis of temperature dependent magnetoresistance, R(7,B), from 
which we deduced: 

1. The charge carriers effective mass, met; 

2. The ground state superconducting coherence length, €(0). 

3. The Fermi temperature, Tr. 

In a result, we find that in all considered scenarios Hx(S,C)y exhibited Te = 190 K (P= 
210 GPa) has the ratio of T./Tr within a range of 0.011 < T./Tr < 0.018. This means that 
H(S,C)y falls in to unconventional superconductors band in the Uemura plot, where heavy- 


fermions, cuprates, pnictides and all near-room-temperature superconductors are located. 


II. Description of the approach 
Detailed description of the approach can be found elsewhere [29]. In short, we use the 
T./Tr ratio to locate the position of the superconductor in the Uemura plot [30]. Tr is 


calculated by: 
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where h = h/2r is the reduced Planck constant, kg is the Boltzmann constant, a = - © 


BT¢ s 
where A(0) is the ground state energy gap. We deduce the ground state coherence length, 
&(0), by the fit of experimental R(7,B) curves to analytical approximate of Werthamer, 


Helfand and Hohenberg theory [31,32] proposed by Baumgartner et al [33]: 


bo po | 


Be2(T) = (2) 


2-1-&2(0) 0.693 


It should be noted that we define T. and B.2 by employing the most strict criterion, i.e. 


R(T) 


> 0 (detailed discussion of the problem can be found elsewhere [34]). We will 


Rnorm 
designate this model as B-WHH model. 
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Because a = cannot be deduced from available experimental data, we perform our 


calculations by assuming that a has lower and upper limits within s-wave superconducting 
gap symmetry: 

3.5 <æ < 4.5 (3) 
where the lower limit is the weak-coupling limit of Bardeen-Cooper-Schrieffer theory [35] 
and the upper limit is the value for anharmonic phonons computed for precursor compound 
H3S by first principle calculations approach [19,36]. 

The charge carriers effective mass, m% ff has been calculated by the use of the 
Eliashberg’s theory expression [37]: 

Mere = (1 + Ac-pn) ' Me (4) 
where me is the electron mass and the electron-phonon coupling constant Ae-pn is deduced by 
the use of the advanced McMillan equations [34,38]: 
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where To is the Debye temperature, and 
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fx = 1 + (0.0241 — 0.0735: u"*) > 2b-pr: (7) 
where u* is the Coulomb pseudopotential parameter (ranging from u¥* = 0.13 0-0.16 
[19,34,38]) for which in this paper we use an average value of u* = 0.13. And the Debye 
temperature was deduced from the fit of R(T, B = 0) data to Bloch-Griineisen (BG) equation 
[39,40]: 
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where To, Ro and A is a free-fitting parameter. It should be noted, that the procedure to deduce 
the electron-phonon coupling constant, Ae-ph, by combined use of the Bloch-Griineisen and 
the McMillan equations is widely used in the field [41,42]. 

It is important to note that analysed R(T,B) datasets for Hx(S,C)y were directly extracted 
from Figs. 1,2 of Ref. 11 (because published R(7,B) datasets are not available from the 
authors of Ref. 11, because these published R(7,B) curves will be in use in a patent 
application [43]). Thus, we extract R(T,B) datasets for the Hx(S,C)y samples by the use of the 
same routine as it has been done for dozens of R(7,B) datasets for other highly-compressed 
superconductors, which we analysed in our recent papers [29,34]. For instance, we can 
mention: black phosphorous (raw data reported by Shirotani et al [44] in their figure 5), 
boron (raw data reported by Eremets et al [45] in their figure 2), germanium arsenide (raw 
data reported by Liu et al [46] in their figure 3), silane (raw data reported by Eremets et al 
[47] in their figure 2), €-phase of O2 (raw data reported by Shimizu et al [48,49,50]), sulphur 
(raw data reported by Shimizu et al [51] in their Fig. 10), lithium (raw data reported by 
Shimizu et al [52] in their Fig. 2), sulphur hydride/deuteride (raw data reported by Einaga et 


al [53] in their Fig. 3(a,b); by Drozdov et al [1] in their figure 2(b); and by Mozaffari et al 


[54] in their figure 1), lanthanum hydride/deuteride (raw data reported by Drozdov et al [55] 


in their figures 1,2,4, and extended data figures 2,3,5). 


III. Results and discussion 

Snider et al [11] in their Fig. 1 reported several R(T,B = 0) curves for three highly- 
compressed H;(S,C)y samples. To be reliably fitted to Eq. 8, the R(T,B = 0) curve should be 
measured at reasonably wide temperature range and also should not to be distorted by any 
experimental artefact (for instance, the failure of electronics or diamond anvil cell, or both), 
which change the shape of R(7,B) curve, and more likely observed drop in resistance, which 
can be mistakenly interpreted as the superconducting transition. We find, that from 9 
presented in Figs. 1,2 [11] experimental R(T,B=0) curves, the only one dataset designated as 
“Run 3” at P = 210 GPa (Fig. 1) can, with satisfaction quality, fit to Eq. 8. It should be noted 
that fits of other datasets to Eq. 8 either do not converge, either after the converging, deduced 
values have some large uncertainty boundaries, exceeding deduced values. This is a very 
unusual behaviour of R(T,B=0) datasets for Hx(S,C)y samples reported by Snider et al [11], 
because 27 R(T,B=0) datasets for a variety of highly-compressed superconductors (ranging 
from elements to superhydrides/superdeuterides) reported by different research groups [1,44- 
55] were fitted to Eq. 8 with a good quality and accuracy [29,34]. In this paper, we analyse 
raw R(T,B) dataset for Hx(S,C)y sample (“Run 3” at P = 210 GPa (Fig. 1) [11]) in an 
assumption that this dataset is accurate and correct. However, it should be stressed that 
independent confirmation/disprove primary results reported by Snider et al [11], including 
very high-T; values, is required [56]. 

In Fig. 1 we show the fit of R(T,B = 0) curve (Run 3, P = 210 GPa) to Eq. 8. Deduced 
Debye temperature is Tg = 1497 + 8 K. By the use of Eqs. 5-7, and taking in account Te = 


190 K, one can calculate Ae-pn = 2.0. This value is in a very good agreement with computed 


Ae-pn = 1.84 (200 GPa) and Ag_ypn = 1.71 (250 GPa) reported by Errea et al [17] for 
precursor binary compound H3S. Based on pivotal result of Duan et al [54] who calculated 
Ae-pn = 2.19 (200 GPa) for H3S, we can conclude that deduced by us Ae-pn = 2.0 for 
Hx(S,C)y (Run 3, P = 210 GPa) is very close to Ae-pn for binary H3S compound. This results 
may allude that sample “Run 3” at P = 210 GPa (Fig. 1) [11] can be just binary H3S 
compound. Based on deduced Ae-pn value, the charge carriers effective mass is: 

Mere = (1 + Ac-pn) * Me = 3.0: me (9) 
It should be noted that this value is in a good agreement with m% sf = 2.76 * Me reported by 


Durajski [16] for the precursor compound H3S. 


zs 
H,(S,C), (P = 210 GPa) 
2.8 
& 
7 2.1 
§ To = 1497 +8 K 
B 1.4 Ae-ph = 2-0 


; ə raw A(T) data 
0.7 ; BG fit 
y @ T; 


150 180 210 240 270 300 
temperature (K) 


Figure 1. Resistance data, R(T), and fit to BG model (Eq. 8) for Hx(S,C)y sample 
compressed at P = 210 GPa (raw data is from Ref. 11, where the sample designated as Run 
3). The fit quality is R = 0.9998. 95% confidence bar is shown. Green balls show bounds for 
which R(T) data was used for the Eq. 8 fit. 

Snider et al [11] in their Fig. 2(b) reported B.2(T) data for this sample. We extract B.2(7) 
data directly from Fig. 2(b) of Ref. 11. The fit of B.2(7) to Eq. 2 is shown in Fig. 2(a). 


Deduced ground state coherence length is €(0) = 2.39 + 0.04 nm. It should be noted that 


Eq. 2 is a good extrapolative tool, that can be proved for the case of highly-compressed 
LaHio. Drozdov et al. [9] reported first B.2(T) dataset for this NRT superconductor for 
applied field up to Bappi = 9 T. This dataset was used to extrapolate B.2(0) value for LaHio in 
Ref. [57]. Recently, Sun et al [58] report new experimental B.2(T) dataset for LaH1o0, which 
was measured on the world-top magnetic field facilities with applied field up to Bappı = 60 T. 
It can be seen that new experimental data is pretty much reproduced extrapolated values 


obtained by the use of Eq. 2 for dataset measured up to Bappi = 9 T in Ref. 9. 
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Figure 2. Superconducting upper critical field, B.2(7), and data fit to B-WWH model (Eq. 2) 
for two Hx(S,C)y samples compressed at P = 210 GPa (raw data is from Ref. 11). The fit 
quality is R = 0.9993. 95% confidence bars are shown. 


From deduced &(0), mý% ff and Te, we calculate value ranges for Tr and Tc/Tr, by 
assuming o within its upper and lower limits (Eq. 3). The results are in Table 1. 


Table I. Deduced and calculated parameters for carbonaceous sulfur hydride. For 
calculations we use deduced me¢¢ = 3.0 * me. 


Pressure Deduced Te Deduced (0) Assumed Tr (10* K) T/T 
(GPa) (K) (nm) 2A(0) 
Kp Te 
1.05 + 0.03 0.018 + 0.001 
210 189.8 + 0.5 2.39 + 0.04 1.74 + 0.06 0.011 + 0.001 


Analysed sample of highly-compressed Hx(S,C)y are shown in the Uemura plot (Fig. 3). It 
can be seen, that in all considered scenarios the carbonaceous sulfur hydride has 0.011 < 
T./Tp < 0.018 and falls in the unconventional superconductors band, where heavy-fermions, 
pnictides, cuprates and highly-compressed near-room-temperature superconductors are 


located. 
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Figure 3. A plot of Te versus Tp where Hx(S,C)y sample (red/black squares) is shown 
together with main superconducting families: elemental superconductors, heavy-fermions, 
pnictides, cuprates, and near-room-temperature superconductors. References on raw data can 
be found in Refs. 29,30,58,60-63. Boundary lines for the Bose-Einstein condensate (BEC), 
the Bardeen-Cooper-Schrieffer (BCS) superconductors, and for ratio of Tc/Tr = 0.05, 0.01 are 
shown for clarity. 


There is also a need to clarify that the London penetration depth, A(T), for which Snider et 
al [11] calculate the value of A(0) = 2 nm, cannot be derived in a way presented in Extended 


Data Fig. 3,b [11]. The estimation for A(0) can be made by taking in account that practically 


all non-elemental superconductors are type-II superconductors, and most of them have the 


A(0) 


Ginsburg-Landau parameter «(0) = FCO) 


within a range of 60 < x(0) < 120 [64-71]. Based on 


deduced €(0) = 2.4 nm, expected value for A(0) is within a range: 144 nm < A(Q) < 290 nm.It 
should be noted that binary highly-compressed H3S has A(O) = 188 nm [71]. 

The amplitude of the ground state superconducting gap, A(0), also cannot be calculated in 
a way, as proposed by Snider et al [11]. The estimate can be obtained by the use of Eq. 3, 
from which one can obtain: 29 meV < A(0) < 37 meV It should be noted that temperature and 


field dependences of superconducting gap, A(T,B), cannot be calculated by the equation: 


A(T) = A(0) - 1.76 + kp «Te: (1 2 =) (10) 


which was used by Snider et al [11] in their Extended Data Fig. 3,d. An accurate analytical 


expression for A(7,B=0) for s-wave symmetry has been given by Gross-Alltag et al [72]: 


A(T) = a(0) tanh | h£- E- = A(O) tan] n€ -1| (11) 


where AC/C is the relative jump in electronic specific heat at Te, and 7 = 2/3. 


IV. Conclusions 

The discovery of superconductivity in highly-compressed H3S by Drozdov et al [1] 
heralded the era of room-temperature superconductivity. To date, several 
superhydride/superdeuteride superconductors have been synthesised: BaHı2 [25], PrHo [73], 
ThHo/ThHio [2], YH4/YH6 [6,7], LaHio/LaDi1 [8,9,14], Hx(S,C)y [11]. The latter compound, 
in accordance with recent report by Snider et al [11], exhibits Te = 280 K at pressure of P = 
267 GPa, which is 5 K higher than T. = 280 K in LaHio (P = 195 GPa) reported by 
Somayazulu et al [8]. In this paper we extract data reported in Figs. 1,2 of Ref. 11 and deduce 
the electron-phonon coupling constant, Ae-ph = 2.0, the ground state coherence length, €(0) = 
2.20 + 0.09 nm, and the Fermi temperature, Tr = (1.7-2.8) 10? K, for the H,(S,C)y 
superconductor with Te = 190 K. Deduced values and reported T. are indistinguishably close 


to values reported previously for primary binary compound of H3S. From deduced values, we 


9 


calculate the T./Tr ratio in Hx(S,C)y compound (T; = 190 K) and find that in all considered 
scenarios this compound should be classified as unconventional superconductors. 

We should stress that our analysis reveals significant problems with all reported by Snider 
et al [11] R(T,B) data for which claimed T: > 200 K. These findings are in a good accord with 
results of independent analysis of R(7,B) data reported by Hirsch and Marsiglio [56]. Based 
on this, independent confirmation/disprove high-T. values in highly-compressed Hx(S,C)y is 


required. 
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